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On 1-loop diagrams in AdS space and the random disorder problem 

Ling-Yan Hung, Yanwen Shan^j 
Perimeter Institute for Theoretical Physics, Waterloo, Ontario N2L 2Y5, Canada 

We study the complex scalar loop corrections to the boundary-boundary gauge two point function 
in pure AdS space in Poincare coordinates, in the presence of a boundary quadratic perturbation 
to the scalar. These perturbations correspond to double trace perturbations in the dual CFT and 
modify the boundary conditions of the bulk scalars in AdS. We find that, in addition to the usual 
UV divergences, the 1-loop calculation suffers from a divergence originating in the limit as the 
loop vertices approach the AdS horizon. We show that this type of divergence is independent of the 
boundary coupling, and making use of which we extract the finite relative variation of the imaginary 
part of the loop via Cutkosky rules as the boundary perturbation varies. Applying our methods 
to compute the effects of a time-dependent impurity to the conductivities using the replica trick 
in AdS/CFT, we find that generally an IR-relevant disorder reduces the conductivity and that in 
the extreme low frequency limit the correction due to the impurities overwhelms the planar CFT 
result even though it is supposedly 1/N 2 suppressed. Comments on the effect of time-independent 
impurity in such a system are presented. 
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I. INTRODUCTION 

AdS/CFT correspondences-OS] has proven itself an extremely powerful tool in extending our understanding of 
strongly coupled quantum field theories, stretching its influence into many different realms of physics. Particularly, 
there has been a surge of interest in applying these techniques in condensed matter systems in recent days (see for 
example a recent review [1] and references therein). Thus far, studies of AdS/CFT have predominantly focused on 
extracting the leading large N physics in the CFT via (semi)-classical supergravity computations in AdS. However, if 
one attempts to make contact with more realistic systems, N should be hnite and there exist many circumstances in 
which l/N correction is important. These corrections correspond to non-planar diagrams in the CFT, and quantum 
loop corrections in the AdS bulk, and in particular 1-loop diagrams in AdS give rise to l/N 2 suppressed corrections. 
A number of physical phenomena are known to show up only if one includes loop corrections in the bulk, such as 
the de Haas - van Alphen quantum oscillations [SJ [6] in strongly coupled charged systems, hydrodynamic long-time 
tails of a fluid [7j , and the holographic manifestation of the Mermi- Wagner therem [5] . More recently it is also shown 
in [S] circumstances where the loop corrections can compete with the tree contribution in certain non-fermi-liquids. 
Consequently, there has been growing interest in the community in understanding quantum loops in the supergravity 
theory in AdS. Calculating loops in AdS space could entail extra complexities, such as additional divergences, due to 
the non-trivial properties of the geometry. It certainly is of importance to investigate this problem closely, clarifying 
possible obstacles and extracting physical implications. 

As one of the original motivations of the current study, it has been proposed in [10] that one may incorporate the 
replica trick that is commonly adopted in the context of condensed matter systems to capture the effect of random 
disorder, through AdS/CFT correspondence. There, one introduces replicas of the AdS spaces and the bulk fields. 
The coupling to the random disorder appears on the AdS side as boundary terms that couple to the replicated bulk 
fields, thus effectively changing the boundary conditions those fields satisfy. These kind of boundary perturbations 
have been studied previously for example in |11H17| . and they are shown to correspond to multiple trace, and in case 
of a quadratic boundary term double trace perturbations, in the dual CFT theory. It is also shown there that the 
effect of these boundary perturbations of a scalar field for example, are only mediated to other fields, such as the 
graviton or the photon, beginning at 1-loop order. This implies that the effect of disorder on transport coefficients 
such as conductivities only shows up in loops 18J. 

Motivated by these interesting studies, we take a phenomenological approach and consider charged scalar fields 
coupled to U(l) gauge fields in pure AdS space in a general Einstein-Maxwell theory in d + 1 dimensions, and study 
particularly the scalar 1-loop Witten diagrams correcting the photon boundary-boundary 2-point function. We then 
read-off the conductivity of the dual theory via the usual AdS/CFT dictionary. These scalars are allowed to satisfy 
mixed boundary conditions, corresponding to the introduction of double trace boundary perturbations as mentioned 
above. It is found that the computation done in Poincare coordinates in Lorentzian signature suffers from a divergence 
that arises as both vertices in the loop approach the horizon (or the IR limit) where the geodesic distance between 
them vanishes. It is surprising, however, that such a divergence persists even within the imaginary part of the loop. 
More specifically, it manifests itself as a singularity in the integral over loop-momenta in the collinear limit j 19], that 
cannot be cured even though the phase-space volume approaches zero there. These divergences, however, turn out 
to be independent of the strength of the boundary perturbation. We demonstrate this property by studying the 
asymptotics of the propagators in the near horizon limit and extract the precise forms of these singularities, even 
though the full analytic result of the loop integral is not generally expressible in terms of elementary functions whose 
properties are otherwise obscure. We obtain in this way both the coefficients and the powers of each singular terms 
analytically and find that they have simple dependence on the number of bulk dimensions. 

Given the universal nature of these divergences, we extract finite results by considering the differences of diagrams 
evaluated at different boundary couplings. We managed to compute the relative conductivities for a wide range of 
boundary couplings / and find that they interpolate smoothly between the two conformal limits / — > and / — > oo. 
We then return to the study of random disorder via the replica trick, and compute the conductivities under the 
influence of random impurities. At d = 2 + 1, we give an example where the computation can be done exactly. We 
find that the presence of the impurities generally reduce the conductivity. Moreover in the low frequency limit the 
correction overwhelms the planar CFT result. This implies that a re-summation in the deep IR is probably necessary, 
but also suggests a possible resolution to the puzzle of how l/N suppressed corrections could actually drastically 
change the IR behavior of the transport coefficients. A more complete analysis of these re-summation is however not 
pursued in the present paper. 

The organization of the paper is as follows. In section III] we present the form of scalar bulk-to-bulk propagators 



satisfying general mixed boundary condition at the AdS boundary. In section III we briefly review the vector bulk- 
to-boundary propagator in momentum space and the scalar-vector vertices that are relevant for the loop diagrams. 
We also introduce Cutkosky rules in AdS space and compute the imaginary part of the scalar 1-loop correction to the 
photon boundary-boundary 2-point function. We did not include similar contribution of other fields such as gravitons 



because we are mainly interested in the leading dependence on the coupling of the boundary perturbation. We will 
show that a divergence arises but is independent of the boundary perturbation, and that the remaining finite part 
interpolates smoothly between the two conformal limits. We also make some comments on these divergences more 
generally in other loops in AdS space. In section |IV| we apply our method in the context of random disorder. We 
conclude our results in section [Vj Further details of the computations are relegated to the appendices. 

II. SCALAR PROPAGATORS OF MIXED BOUNDARY CONDITIONS 

Consider AdS space in d + 1 dimensions with the metric 

ds 2 = -^{-dt 2 +dz 2 +^2 dx 2 ), (1) 



the Green's function of scalar fields in AdS space of mass m satisfy the inhomogenous Klein-Gordon equation sourced 
by delta-functions: 

(-^(d Zl g zz Vgd Zl )+g^d^-mAG A _(z 1 ,z 2 ,x 1 ,y 1 ) = - [ -5 d (x 1 -y 1 )S(z 1 -z 2 ). (2) 

In general these scalar fields behave in the boundary limit z — > as 

cf>{z) ~ az A + + f3z A - , A± = ^±v, (3) 

where we have defined 



d 2 
v=\jm 2 + — . (4) 

When both A± are greater than zero and thus the corresponding wave-function normalizable, one can insert boundary 
terms of the form 

SS aM = [ f(3 2 , (5) 

JdM 

so that the variation of the action only vanishes on the boundary if the scalar field satisfies the following boundary 
condition 

a = f/3. (6) 

For finite / the dual operator in the CFT whose vev is given by /3 has conformal dimension A_ in the UV, and 
therefore the dimension of / is given by 

[/]=d-2A_=2i/. (7) 

The unitarity condition and the requirement of preserving conformal symmetry in the UV, i.e. the perturbation is 
irrelevant in the UV, demand \v\ < d/2, or -d 2 /4 < m 2 < [217] . 

These boundary terms are multi-trace perturbations from the point of view of the dual CFT, and have been first 
studied in e.g. pdUT7] . 

The scalar Klein-Gordon equation can be readily solved in momentum space of the flat directions. For given d- 
momentum k^, x ' 2 J v (kz) is a solution to the homogeneous version of equation (pi) i.e. without the (5-function source, 
where J u (kz) is a Bessel function of order v and k = |fc p |. One can easily solve for the Green's function making use of 
the completeness of Bessel functions for given boundary condition at z — > 0. Instead of piecing two solutions defined 
for Z\ > z 2 and Z\ < z 2 as in [21] , we can express Ga_ ( z ii z 2, k) as the following integral 22 

f°° A 

GA_(zi,z 2 ,k) = dA J v j(A,zi)J u j(A,Z2), (8) 

where 

J vJ (A,z) = N(A)[A(A)J v {Az) + B(A)J_ v (Az)] (9) 



is the linear combination of both Bessel functions of order u, properly chosen to satisfy the boundary conditions ^ 
at z — >• 0. It is straightforward to check that we must choose 

a = i, B ^^l^A v ( 10 ) 

and the corresponding normalization factor 

N ^ = 1 + 2B{A)cos(vtt) + B(A) 2 ' ^ 

such that J v j is normalized in the sense that 

zA 1 J J/if {A 1 ,z)J v j(A 2 , z)dz = S(Ai - A 2 ) . (12) 



o 



The expression for N(A) can be justified most easily by checking the asymptotic forms of J v and ./_„, or using the 
orthogonality properties of the standard Bessel functions. 

One could check, as we demonstrate in AppendixIXJ that this representation agrees with the bulk-to-bulk propagator 
in [21) obtained in the Euclidean signature, which simply means k 2 is positive in J8J). The virtue of making use of this 
representation will be made manifest when we begin computing loop corrections, where we are spared of the difficulty 
of dealing with the step function 0(z — z') present in the representation in [21] . Since we are interested in Lorentzian 
signature we will have to specify precisely the causal structure of our propagators. In the following analysis, we will 
consider loop corrections to a retarded correlation function. Feynman, or time-ordered, propagators, as discussed in 
|23j . can be obtained by a simple is prescription. Namely, one makes the replacement 

(13) 



A 2 + fc 2 A 2 + k 2 -ie' 

Similarly when we consider retarded (advanced) propagators, we will put all the poles of fco in the lower (upper) half 
complex fco plane. 

III. 1-LOOP CORRECTION TO GAUGE TWO-POINT CORRELATION AND CUTKOSKY RULES IN 

ADS SPACE 

We are interested in this paper the charged scalar 1-loop correction to the boundary-boundary correlator of a 
U(l) gauge field in AdSd+i. According to the standard AdS/CFT correspondence, a 17(1) conserved current in a 
d-dimensional CFT is generally dual to a U(l) gauge field, or the bulk "photon" in a d + 1 dimensional AdS space. 
One can view that the CFT lives on the AdS boundary and there is the relation 

We will be considering, on the AdS side, loop corrections to this correlator from a charged minimally coupled scalar 
field that satisfies the general mixed boundary condition explained above. On the CFT side, these loop corrections 
correspond to 1/A^ 2 correction to (jfj,j v )- 

A. Interaction vertices and the photon propagator 

Let us first briefly review the photon boundary-to-bulk propagator at tree level and its interaction with the complex 
scalar. For convenience, we will work in the gauge 

A z = d^A" = 0. (15) 

In Euclidean signature, the boundary-to-bulk propagator satisfying the usual Dirichlet boundary condition at z = 
in this gauge is [23H26] 

A (zv)- j± (v] ^ d/2 ~ lK W2-i)(P^ „ J±() _ Q ri „v 

A,(z,p) - J^ (p) {pe)d /2-i K{d/2 _ i){pe) > P J , IP) °- ( 16 ) 



The propagator has been normalised and e is the UV (or boundary z — \ 0) cut-off. The boundary source J- 1 satisfies 
the transverse condition as a result of the gauge choice. The Lorentzian propagator can be readily obtained from the 
Euclidean one by analytic continuation. Depending on whether one is interested in in-coming or out-going boundary 
conditions at the horizon z — > oo, one could replace p by ±ip accordingly. Since we are computing a retarded 
correlation function, we will take one external leg to satisfy in-going boundary conditions, whereas the other leg 
should be the complex conjugate [27] . 

There are two types of interaction vertices between minimally coupled charged scalars and photons. The 3-point 
and the 4-point vertices are given by 



V MH = / d d xdz^gg^(-i)A^d v <f>-<f>d^) 



-gg^A.A^ 



Vaa^h = / d d xdz^ 



(17) 



Given that the sources J 1 - in our gauge have to satisfy the transverse condition, only the component satisfying 
the Ward-identities p^ ( j^jv ) = would contribute. For general external momenta p the 1-loop contribution is 
not expected to satisfy the d-dimensional flat-space Ward identities for the boundary theory. When p approaches 
zero however, the photon wave function becomes independent of the radial coordinate and the integral over the 
radial position of the vertices yield simple delta functions. In this limit the amplitude is reduced to a d-dimensional 
calculation, and can be shown to be transverse after all 1-loop diagrams are combined as usual in flat-space. 

B. 1-loop diagrams and the Cutkosky rules in AdS space 

Now we turn to the 1-loop correction to the boundary-boundary correlator of the bulk photon. Our main objective 
is to find their imaginary contributions. There are two loop diagrams relevant to the current study as shown in figure 
(fTl). Both diagrams contribute to the 1/N 2 correction to the current-current correlation on the CFT side. Let us 





(a) 



(b) 



FIG. 1. 



denote the contribution from the diagram (a) as {j^j„ ) . Using the propagator (8) and (63), this diagram boils down 
to the following integral 



(M-P)3u(p))' = £d _ 2 | A , 



P m (p)Pvp{p) f d'fcdAidAa (p- 2k) a (p- 2% A 1 A 2 \I vJ (ip,A 1) A 2 )\ i 



d/2-l 



(*pc)| s 



(27T) 



(A? + F)4AI + (p _ fc)2 )fl 



(18) 



where we have included the only non- vanishing contribution to this loop, involving the product of a retarded and 
an advanced scalar propagator. Our readers are reminded of the respective pole structure of the propagators by the 
subscripts R and A above. We imposed the gauge condition by introducing the usual projection operator 



We have also defined the "vertex function" I v j as 



Vp,Vu 



poo 

I u j(j3,Ai,A 2 ) = / dzz^Ka_ 1 (pz)J l/ j(A 1 ,z)J„.f(A 2 ,z), 
Jo 2 



(19) 



(20) 



which can be read-off directly from equation (17) expressed in momentum space for the d flat directions. 

More generally, it is possible that the two scalars propagating in the loop satisfy different mixed boundary conditions, 
in which case one would have to generalize the "vertex function" in the loop-integral as 

/>oo 

I",fi,h(P:M,M)= dzz^Kd_ 1 (pz)J Vl f 1 (A 1 ,z)J v j 3 (A 2 ,z). (21) 

Jo 2 

This function has an analytic expression that will be discussed further below. At the moment, we just note that, 
following immediately from the properties of the Bessel functions, 



L 



■ vJl , h {e^ P ,A 1 ,A 2 )=Il Juh {e-^p,Ai,A2), (22) 



which we have already made use of in equation ( 18 1 



Clearly, the integral would suffer from divergence near the boundary z — > 0, interpreted as a UV divergence, if 
2\i/\ > 1. To avoid complications, we will restrict our attention in this paper to \v\ < |. 



Written in this form, integral ( 18 ) can be conveniently interpreted in the d-dimensional language as the total 
contribution to the 1-loop correction of ( A^A^ ) from all pairs of KK-scalars of masses Ai and A2. These KK-scalars 
form a continuous infinite tower whose mass take all values from to infinity, as should be expected since the z- 
dimension is non-compact. Their contribution to the loop is integrated over, convoluted by a non-trivial "vertex 
function", given by I v j(ip, Ai, A2), or more generally Iv,fi,fa describing the mixing of the KK-scalars and their 
coupling to the gauge field. It should be mentioned that the vertex function I v j 1 j 2 as defined in (20) is not a proper 



integral and needs to be regularized. Most naturally, it is regularized by adding to the external momentum p a small 



positive imaginary part ie so that (20) becomes well-defined since its integrand decays exponentially as z — > +00, and 
then taking the e — > limit. The result is in general finite and well-defined, but may contain various singularities 
depending on the values of (p, A 1; A 2 ). Some of the singularities are severe enough that they can render the loop 



integral ( 18 1 ill-defined. Below, we will devote a significant part of this paper discussing these singularities and the 



regularization of the loop integral. 

In this d-dimensional language, one can readily obtain the imaginary part of the loop ( 18 ) using the Cutkosky 
Rules, or the analog of the "optical theorem" for S-matrices. At the level of tree-diagrams this has been considered in 
[23] . Analyzing the pole structure for the product of retarded and advanced propagators in the same vain shows that 
the imaginary part can be evaluated by putting the propagators in the loop on-shell, and, as in Feynman loops, only 
the physical poles would contribute. The cut diagram corresponds to the amplitude of a single "photon" decaying into 
two KK-scalars of masses Ai and A2 respectively. One then finds the modulus of the amplitude and integrates it over 
the phase space allowed by the kinematics as well as all positive values of Ai and A 2 . Equivalently, one can replace 
the scalar propagators [(Ai+k 2 )(A 2 + (p — k) 2 )]^ 1 by their on-shell conditions, i.e. — (27r) 2 5(A 2 + k 2 )S(A2 + (p— k) 2 ), 
and carry out the remaining integral. The presence of the two <5-functions confines the integral into a compact 
region in the phase space. In fact the result is only non-vanishing if the incoming momentum p is time-like and 
x /pg-|p| 2 >A 1 + A 2 . 

Consider the simple example where p^ = (po, 0, ..., 0) and evaluate the spatial components of the loop diagram. By 
rotational symmetry and via the Cutkosky rules, we find 

pSb *<M(a! + (p- kf)S{k\ + (if) 

8(d- 1)(4tt) c 



Sij 



*« «m _ 1 a° w-2 ft - (Ai + A 2 ) 2 )(l - (A 1 - h) 2 ) , (23) 



where fl m is the volume of an m-sphere and A1.2 = Ai^/p. The time-time component of the momentum integral 
similarly evaluates to 

j^^iPo 2k a ) 2 S(A 2 2 + (p- k) 2 )5{A\ + (k) 2 ) 

n d -2P d - 2 ,72 A^2/n ,1 , 1 ,2.1 ,1 a ^2^^ (d - 3)/2 



8(4 



^(A?-A^((1-(A 1 + A 2 ) 2 )(1-(A 1 -A 2 ) 2 )) " . (24) 



Putting the pieces together, we are left with the integrals of Ai and A 2 only. By the AdS/CFT dictionary, we find 

Sjj n d _ 2P d + 2 

8(d- l)(47r) d - 2 e^2|A> 2 _ 1 ( l p e)| 2 



T / • I \ ■ / \ \' $ij Qd-2P 



where we defined 



( 25 ) 

■ / dA x / 1 dA 2 A 1 A 2 iJ(A 1 ,A 2 )( d - 1 )/ 2 |^ / ( i , A X ,A 2 )| 2 , 
Jo Jo 



H(x,y)= [l-{x + v) 2 } [l-(x-y) 2 ] (26) 



and will refer to it as the "phase volume" factor. The correlator ( jo(— Po)jo(po) ), on the other hand, is killed off by 
the projector -P M „(p) for the particular external momentum (p , 0, . . . , 0) we have chosen here. 

As a consistency check, consider the simple case where the boundary term is set to zero i.e. / = 0, where conformal 
symmetry is expected to be preserved. In that case the po dependence of I v j can be completely taken out, which is 
simply 1^^ ~ p~ 1_d ' 2 . Using also the fact that lim e _>o Kd_ l (ep) ~ {ep)~ 2" 1 , the correlator is therefore given by, 

(M-Po)3i(po)Y = Si3&t 2 , (27) 

where £ is a p independent constant obtained from the rest of the integral. The po dependence of the correlator is 
exactly what is expected from the conformal dimension of a vector current in a CFT in (^-dimensions. 

In the case where / is non-zero, the simple scaling behaviour is disturbed and conformal invariance broken. It is 
manifest that in this case the external momentum p appears as an additional parameter in the loop integral. 

We make a comment regarding the loop diagram (b) in figure (IT]). It is necessary to include this loop to recover 
gauge invariance in the bulk. But the external momentum in this loop factors out completely and it does not carry 
any imaginary part, as in the case in flat-space. To extract the imaginary part this diagram can be safely ignored, 
and its real part can be inferred via gauge invariance from diagrm (a) . For this reason, we will omit this loop diagram 
in what follows. 

The picture outlined above is mostly standard. One is easily tempted to assume that the loop integral can be 
evaluated without any obstacles. There, however, remains a quite serious difficulty. As we have already alluded to, it 
turns out that the "vertex function" I v .f 1 j 2 is not a regular function for all values of Ai and A 2 . It contains various 
divergences of different orders. Expressed in the 2-dimcnsional space spanned by {A^ A 2 } the region where Ii>j 1 .f 2 
becomes singular is not isolated but form boundaries, and it diverges sufficiently fast near those boundaries that it 
can cause the loop integral to diverge, even after the vanishing of the "phase volume" factor _ff( d_1 )/ 2 in this limit is 
taken into account. This divergences is a new type of divergence, different from the usual UV divergence which still 
comes about from the remaining d k integral. To compute the imaginary part of the 1-loop correction and obtain a 
physical answer, one must first understand the singularity structures of I v ,f 1 j 2 an d regularize the loop accordingly. 
This is the topic we turn to in the next subsection. 

C. The divergence of the vertex function and the 1-loop integral 

Some discussions on the divergence properties of the "vertex function" Iu,f lt f 2 are i n order. We have mentioned 
that as / = 0, 00, the p dependence of I Vi q or 7„ i00 can be completely taken out and apart from an overall power-law 
dependence of p, the rest, including any possible poles of course, depends on the ratio Ai. 2 = Ai^/p only. When 
/1.2 7^ 0, there is no longer such a simple scaling behavior of Ivj 1 .f 2 - However, as it's explained later in this subsection 
as well as in Appendix lB| as far as the divergence of Iuj ± j 2 is concerned, a simple scaling behavior still arises and its 
singularity structures depend on the ratios Ai j2 and / = f/p only. Therefore, without loss of generality, we will set 
p = 1 and assume A 12 = A 12 in this subsection, ignoring the scaling of / whenever it's unimportant. 

We will outline the main point and provide the essential results in this section, and leave the full details in Appendix 
[5] for those who are interested. 

Impressively, Bailey [2 8) evaluated the following integral involving a power and three Bessel's functions and found 
an analytic answer: 

A-i ts ,„\ r ,__w n..-\ - 2 x - 2 a^T[l/2(X + y. + v + p)]T[l/2(X + fx + v - p)\ 



dzz K p (cz)J lJi (az)J v {bz) 



i x +^+"T[y + l]T[v+l] 



Fi[^ + H + v - p); -(A + n + v + p); fi + 1; v + 1; -a 2 /c 2 ; -b 2 /c 2 } 



Here, F4 is the Appell hypergeometric function, one of the generalized hypergeometric function that, apart from four 
parameters fixed in the current case by the coefficients A, p, \x and is, depends on two complex variables. For the 
loop integral that we wish to compute in the Lorentzian signature, we need to analytically continue Bailey's result, 
substituting c by ic, and obtain the following identity: 



I 



+ OG 



z 2 K d_ 1 (iz)J tl (kiz)J v {k 2 z)Az 
T[l + \{ti + v)]T[\{d + p + v)\k>tk» 2 

(i)«*/2+i+^+«T(/i + i)r> + 1) 



(28) 



Fi[l + ~(M + v), -(d + \x + v); n + l,v + 1; A?, k\ 



This equation, obtained by analytically continuing the convergent integral ( 28 1 , agrees with the regularization scheme 
for Iv,fi,f 2 explained earlier where one endows p with a small positive imaginary part that is taken to zero in the end. 
Within the region |Ai| + |A 2 | < 1, the Appell hypergeometric function has a series expansion which fails to converge 
as |Ai| + IA2I — > 1. A simple observation of the asymptotics of the integrand in (28) when z — > 00 leads to the same 



conditions for / to be well defined. In the current context, both Ai^ are real and positive, which allows us to consider 
them as the mass of the KK-scalars, and the condition Ai + A2 < 1 appears to be nothing other than the consequence 
of energy- momentum conservation, i.e., a particle can only "decay" into two scalars whose masses add up to a value 
smaller than the proper energy of the decaying particle. But the fact that / approaches infinity as Ai + A2 — > 1 is 
somewhat problematic and a regularization scheme is needed. 

To this end, we need to understand exactly how the Appell hypergeometric function F4 diverges as A1.2 approaches 
the convergence boundary. It is sufficient to investigate the asymptotics of the integrand in ( 28 1 as far as the 
singularities are concerned. We assume that the values of d, v and fi are chosen properly such that the integrand is 
regular at any finite value of z, or if it does contain any singularity at z < +00, it is a singularity that can be integrated 
across and leads to no singular behavior of I. Given such assumptions, the only possible source for / to diverge is 
when the integral is carried out all the way toward z — > +00. In the region that z is large, we can approximate the 
Bessel functions by their asymptotic expansions and formally write 



1 



00 n n—K 

EEE 



- ^-l,n-k-l){p,k){v,l) 
dz — ^ 

k \l 



2 5 ^A7A2^ ^ t^° (2zr+^AfA' 2 

id the chips: 
form is given in (B7). The notation (y, n) is defined as 



f e -i(6n+6»+ 3 l ) i 2(k+l)-n e -i[l-(A 1 +A 2 )]z + _ _ _ \ 



(29) 
where 9^ = ^ and the ellipsis represents similar terms that involve other combinations of 1 ± Ai ± A2, whose full 



(u,n) 



nira + is-n) 



(30) 



If we just carry out the integration term by term, regularizing the oscillatory integrand by lifting the integral contour 
slightly above the real axis as explained earlier, we obtain: 



I 



y^ y^ ^ fc (f-l,n- k - l) (/,, k)(u, l)T( 



d-1 
2 



2™AfA 2 



(31) 



^ e -i(flM+^+^f^)(_)fc+»[i _ ( Al + A 2 )] n+ ^ + . . . ) 



Again, the ellipsis represents similar terms that involve other combinations of 1 ± Ai ± A2 and are given in full in 
|B9| and |BTo|. 



We must elaborate a bit more on the procedure outlined above. The two formulae just given would be flawed if 
the "~" were taken to be equal, because the expansions on the r.h.s. of these "equations" are asymptotic expansions 
only. At any fixed value of z, the infinite sum does not converge. Usually, integrating an asymptotic expansion term 
by term is only meaningful if the result is also considered as the asymptotic expansion of the true integral when the 
lower integral limit approaches infinity. However, since for a given integer N, the sum of the first N terms in the 
expansion approximates the full integrand with an arbitrarily small error when z — > 00, if one formally integrates 
the expansion term by term and drops everything that remain finite as |Ai| + | A.2 1 — > 1 along any trajectory, only 
a finite number of singular terms remain and they describe precisely the same singularities of the original integral 
(28). Readers may find it suspicious that integrating the asymptotic expansion from z — to +00 could lead to 



anything meaningful since it is only a good approximation to the actual integrand when z is sufficiently large. More 
appropriately, one should choose a cutoff scale L that is large but fixed and separate the infinite integral into two 



9 



parts: an integral from to L, and an integral from L to +00. The first piece necessarily contains no poles of Ai i2 
because L is finite and the integrand is regular. The second piece must consequently include all the singularities of 
(28), which are independent of the arbitrarily chosen cutoff scale L. [29. Therefore, as far as the singularities are 



concerned, one can choose an arbitrary cutoff L > 0, replace the integrand by its asymptotic form, and carry out the 
integration from z = L to 00 term by term. The resulting expansion contains finite number of terms that diverge as 
I Ai| + IA2I — > 1. Evaluated sufficiently close to the poles, these divergent terms become independent of L, and setting 
L = is only the most convenient choice. We must emphasize that this method would not in general lead to any 
useful information regarding the finite part of the original integral. Therefore, the "~" sign in equation (29) means 
"equal up to an arbitrary regular function" . 
It's worth noting that exceptional cases do exist when such an analysis leads to stronger results. In particular, if 



in (28) d is odd and /i and v are half integers, the asymptotic expansions for Bessel functions of half-integer orders 
are truncated to finite sums in which case the "~" can be replaced by the equal sign and the equations are exact. 
Consequently, for d = 3, one finds examples where a full analytic result can be obtained, as presented in section IV B 
In the conformal limit, we take [i — v and readily find the singular terms of /„ j0 o or I v $ from equation (31) when 
1 — (Ai + A 2 ) — > 0. For general "mixed" boundary conditions we consider, / ^ and the dependence of I v f x t 2 on 
p is more complicated. But the asymptotic behavior of J„./ as z — ¥ +00 is quite simple and essentially identical to 
that of standard Bessel functions J±„, except for a simple phase shift explained in details in Appendix [Bj We can 
therefore follow the same procedure and easily obtain the most important result of this section: the singularities of 
\Ivj t j 2 \ 2 within the domain < Ai, A 2 < 1 are described by 



l i »',/ 1 ,A(*> A l. A 2)| S 



1 



8^AiA 2 Ul-(Ai+ A 2)] d - 



1 



2(M) 
2 rf-3 

fOf[l-(Ai 



1 

aT + 

A 2 )r 



1 

A^ 

(d-l)/2 



1 



[1-(A!+A 2 )] 
{d > 3) . 



d-2 



(32) 



We've arrived at a pleasant surprise, which is probably physically well expected: the leading first and the second 
divergence of \Iu,f lt f 2 (h A i> A 2)l are independent of the parameter / that describes the mixed boundary condition 
for the scalar field <f>. The leading divergence is furthermore independent of its bulk mass parameterized by v. The 
remaining singular terms, on the other hand, do in general depend both on / and v. But fortunately, the orders of 
those "true" singularities are equal or less than (d — l)/2, just sufficiently low to be suppressed by the phase volume 
factor 



H {d-i)/2 = {[1 _ (Ai + Aa)] [! + Al + a 2 ] [! _ (Al _ Aa )] [! _ ( Aa _ Al )]} (d - 1)/2 , 



(33) 



which contains precisely all four combinations of (1 ± Ai ± A 2 ) to the power of (d — l)/2. Hence, the /-dependent 
imaginary part of the 1-loop integral (25 1 is finite and can be evaluated free of any pathologies. This conclusion holds 
true for any d > 2. 

It should be noted once again that, in (31), only terms that become singular as 1 — (Ai + A 2 ) — > should be kept 
and the rest must be ignored, unless the asymptotic expansion is known to converge. For example, when d < 4, only 
the leading divergent term in (31 ) is really there. The 1/[1 — (Ai + A 2 )]( d ~ 3 )/ 2 term and therefore the second term in 
expansion ( [32] ) do not exist. Similarly, terms as 1/(1 — Ai — A 2 )( d ~ 5 )/ 2 in (31) is to be discarded if d < 6. Precisely 
as the method fails, the coefficients in expansion (32) become singular. 

We remind the readers that for general values ofTy, Iv.f 1 ,f 2 contains other singularities whenever 1± Ai ± A 2 — > 
and an analogous expansion can be obtained in the same manner near each of them as shown in Appendix [B] [SU\- I n 
fact there's a natural way to understand why all those singularities are essentially the same. We take this chance to 
mention some other interesting properties of the vertex function which might be useful for a full calculation of 1-loop 
integrals in AdS space. For brevity, we will ignore the differences between photon and scalar and consider a bulk to 
bulk 1-loop integral [3T], in which the following vertex function appears: 



/(P.A] 



■ A 2 )= / 
Jo 



dz z a J v (pz)J u (k 1 z)J l ,(K 2 z) . 



(34) 



where a and v depend on the details of the theory. We'd like to point out that if one can ignore the subtleties related 
to the singular behavior of the integrand at z — > 00, formally this function / is defined on CP 2 if all variables are 
complex numbers, apart from a trivial power-law factor. This is because one can always absorb an overall scaling of 
(p, Ai, A 2 ) — > (Xp, AAi, AA 2 ) by redefining z — > \~ 1 z together with a change of integral contour which always connects 

to 00 in the complex plane. If we restrict ourselves to the real space, one can only absorb such a scaling if A > 0, so 

1 is defined on S 2 . Let us focus on the case when p and Ai 2 , are real and positive. A simple scaling of z immediately 
allows us to set p = 1. We find there are four different regions in the first quadrants in the (Ai, A 2 ) plane separated 
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from each other by singular boundaries. They are I = {0 < Ai + A2 < 1}, II = {1 < Ai — A2}, III = {1 < A2 — Ai}, 
and IV = {1 < Ai + A2 , |Ai — A2I < 1}- We have been focusing on region I only in the above discussion. But region 
II and III can be mapped to I by a simple scaling. To go from region II to region I, we just scale z — > AjZ in (34 1 and 
define A 2 = A 2 /Ai, and Ai = 1/Ai. It's easily verified that < Ai + A 2 < 1, and 



/(1,A 1; A 2 ) 



A«+ 1 /(1,A 1 ,A 2 ) 



(35) 



Similarly, one can go from region III to region I by a simple scaling of z — > A 2 z, defining A2 = I/A2 and Ai = A1/A2, 
and verifying that 7(1,A 1 ,A 2 ) = A^ +1 1 (1 , Ai , A 2 ) . Notice that I(1,A 1 ,A 2 ) actually vanishes quite fast as A 12 is 
large. Region IV is more complicated, but if one is only interested in the asymptotics of the integrand in ( 34 ) , which 
dictates the singularities of I as we argued above, one can also map it to region I by a scaling z — > (Ai + A^Jz, define 
Ai = (1+Ai - A 2 )/[2(Ai +A 2 )] and A 2 = (1 + A 2 - Ai)/[2(Ai + A 2 )], and find that the singularities of J(l,Ai,A 2 ) 
are the same as those of 1(1, Ai, A 2 ) but with (Ai, A 2 ) located in region I now. Therefore all the singularities of / as 
I Ai| + |A 2 | — >• 1 are related to each other in a simple way. 

Just as a demonstration, we show in figure (I2| the numerical result for the relative variation of the imaginary part 
of the 1-loop diagram against the boundary coupling / when d = 4 and v = ~ . The numerical result presented in the 
plot is subtracted against that when / = 00. 
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FIG. 2. The relative variation of the imaginary part of the 1-loop integral against coupling / at d — 4, v — |. External 
momentum being an overall scale is set to one. 



D. A few more general observations on the divergences in AdS loops 



Thus far we have studied closely loops involving charged scalars coupled to photons and the divergences associated 
to the vertex. We would like to make a few observations about loops in AdS more generally. We found that this kind 
of divergence occurring as the vertices approach the AdS horizon is fairly general. Consider for simplicity a <p n vertex 
of the form 



. n „oo / n \ n 



(36) 



expressed in momentum space for the flat directions. Substituting in the integral the radial wave-function of the scalar 
field of mass m, which, for Neumann boundary condition as an example, is given by z d ' 2 J^(pz), where p — \p^\, and 
taking the near horizon limit, we find 



/ n / n \ />oo 



COS 



£,=±1 



/XsjPjZ + Sj6 v ) 



(37) 
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We have used the asymptotic behavior lim z _ i . 00 J v (pz) ~ -4= cos(pz + 6 V ). By naive power counting, ignoring the 
oscillatory factors, this integral is only convergent if 

n -^<d + l, (38) 

precisely identical to the renormalizability condition one would find via naive counting of engineering dimensions. 
Notice this analysis is independent of the boundary condition for the scalar since it would only affect the oscillatory 
factors and leave the power-law dependence intact. 

Such a coincidence no longer exists as soon as we consider coupling to fields with higher spins. In the case of 
photons, the vertex would appear divergent by the above naive power-counting for arbitrary d independently of the 
renormalizable conditions, because the coupling invariably involves the inverse metric which provides extra factors 
of z 2 even though the photon propagator is proportional only to z~5~ J d / 2 -i, and subsequently contribute to the 
divergence in the horizon limit. This consideration is generally applicable to the coupling of scalars and graviton as 
well. 

It is very important to note that even as the z-integral appears divergent according to the above consideration, it 
does not, however, necessarily imply that the loop integral is divergent. The vertex integral can be regularized by 
analytic continuation as we have done in our earlier discussions making use of its oscillatory factors. Using the method 
introduced in the previous section, one can easily verify that this vertex after being regularized must include poles at 
St £ iPi = Oj wnere e i — il) an d could potentially lead to similar divergences in loop integrals in the collinear limit 
as Si Pi ~~ ^ 0- Depending on the order of these poles, which must be analyzed case by case following the procedure 
explained in section |IHC| this may or may not lead to pathologies. Let us mention, as an example here, the same 
analysis implies that this extra divergence occurring as the vertex in diagram (b) in figure (nj) approaches the horizon 
is logarithmic if d = 4. 

We would also like to point out here that using the form of fermion wave- function in AdS as given in [35] , we can 
repeat the above power counting for a fcrmion-photon vertex as z — > oo and find 

/OO /"OO 

dz^/^gA^e^T' 1 !/) ~ / dzz^~ x oscillatory factors , (39) 

which is exactly the same as in the case of scalars. Here e^ is a simple diagonal choice for the vierbeins, and T a are 
flat space gamma matrices. It seems possible that in supersymmetric theories these type of divergences are cancelled 
automatically among fermions and bosons. 

Finally, we make a passing comment about potential UV divergence near the boundary z — ¥ 0. We have restricted 

our attention to scalar fields of mass m such that v = U 4- + m 2 < g because this is a window of consistent masses such 
that UV divergence near the boundary is simply non-existent, and curiously this condition is dimensional independent. 
It would be interesting to have a better understanding of the physics behind these observations. 

IV. RANDOM DISORDER AND THE REPLICA TRICK IN ADS/CFT 

We now turn to an interesting application of the results given in the previous section and calculate the effect of 
random disorder on the transport coefficients using AdS/CFT correspondence. 

A. The replica trick and the conductivities 

We begin with a brief review of the replica trick. Consider a quantum system described by the action S. The effects 
of random disorder can be captured by introducing into the action an extra scattering potential SS 

5S= f d d xV(x,t)0(x,t), (40) 

where V(x, t) is a general spacetime dependent random scattering potential, and O is some physical operator, such as 
the charge density etc. In the presence of the random scattering potential term, correlation functions averaged over 
the random variable V is given by 

J D<p e -f^O (xi)0(x 2 ). 

j D<f>e- 



(0( Xl )0(x 2 )...) = J D[V)P[V) ( J " r „^s -ss ) > ( 41 ) 
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where P[V] is the probability distribution of the random potential and the over-line denotes averaging over V. A 
well-known trick in the condensed matter literature to compute this averaged correlator is called the replica trick. 
The idea is to introduce n copies of the theory concerned such that all operators in the theory are replicated n times 
and we label them by an extra index i G {1, ..., n}. The partition function of the full theory including the n copies is 
related to that of the original one by Z„ — (Zi) n . Then the correlation function for a fixed V in the original theory 
is formally related to the n-replica partition function, treating n as a continuous variable by analytic continuation as 



(OO. 



n->o aJi a,)\ n 



J i= Q 



v S 5 l ? 
Iim — -— -...-Z r , 

n-»0 oJ\ oJ\ n 



(42) 



,7;=0 



where Jj is the current coupled to the replicated operator Oj. The averaged correlation functions can be readily 
evaluated as 



n-^-o dJi dJi n J 



Ji=0 



r S S l 9 



Ji=0 



We have exchanged the order of taking the n — > limit and the integral of V in the first step and defined 

-T,i = i(Si+SS,)+lnP[V] 



Z n = f D\y]Y[D[<f>i]e-ZZ- 



(43) 



(44) 



in the second. For more complete explanation, see [33] and references therein. 



simply Gaussian and local in space and time i.e. 



This relation (43 1 is useful because in the simple but common situation where the random distribution P[V] is 

(45) 



V""(<M) 

P[V] = e v- 



the scattering potential V can be readily integrated out, and we have 



n*[ 



,-T^iSt-i j ^^0*0' 



(46) 



To study the effects of random disorder in a strongly coupled theory, it is natural to incorporate the replica trick 
in the context of AdS/CFT correspondence, which is considered in [101 [34]. The spirit of the two studies are very 
different and in this section we follow [TU] and attempt to study the effects of disorder on the conductivity of a general 
charged system. 

In |10j the replicas are introduced literally as n copies of the AdS backgrounds. The random scattering potential 
in the CFT in the replicated theory is presented by the boundary term 



which gives 



SS = 



SS 



0M d 



d d xV(x,t) 






O' 



d' 1 . 



dM d + 



;^20 i O j 



(47) 



(48) 



i -j 



after the potential is integrated out. These boundary terms are precisely those multi-trace perturbations in the dual 
CFT theory discussed in earlier sections, and dictate the boundary conditions of the corresponding scalar fields in the 
multiple AdS bulks. They also relate the boundary values of the scalars in different AdS copies and thus connect them 
together. The modified boundary-boundary and bulk-to-bulk propagators of these scalar fields have been studied in 
[2"T] and, in the context of the replica trick, in [ID] . It is known that the effect of these boundary perturbations 
are mediated to other sectors beginning only at 1-loop order [TJ] [351 I3S] • Using the folding trick that treats these 
replicated bosons as if they live in the same AdS space [3S] we can thus apply techniques discussed in the previous 
sections to compute the conductivities in the presence of random disorder. 

Now consider again scalar field <j> with mass m. Recall that its boundary expansion is given by 



az 



^ 



± 



(P 

T' 
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In the presence of the boundary deformations (48 ), the scalar fields concerned satisfy nontrivial boundary conditions. 

n 

<* = /!>' (49) 

3 

where n is the number of replicas introduced. These conditions can be diagonalized by: 

n _. 

^=$>i&' 0</<n, a " = ^> E a f" = - ( 5 °) 

3 3 

There is only one linear combination of the fields fa (corresponding to <j) n defined above) satisfies the mixed boundary 
conditions depending on the strength of the disorder /, i.e., 

a n = nfJ3 n , (51) 

whereas the rest of the fields fa^ n simply satisfy usual Neumann boundary conditions 

a mi t n = 0. (52) 

It would be useful to pick an orthonormal basis for these fields, so that the bulk-to-bulk propagators take a simple 
form. One convenient choice is 

1 r 1, l<»<(n-0 

V(n-l){n-l + l) I o, i>n-l + l. 

The bulk-to-bulk propagators of the fields fa in the original basis would be related to these rotated basis ij> by 
simple linear combinations. In particular, we are interested in the following propagators 

G nn = -(G nn + (n- l)^"- 1 ^"- 15 ), G' 3 = (fafa), G 11 = (^), 

Qn(n— 1) _ _(Qnn _ Q(n—X){n—l)\ (54:) 

n 

The bulk-to-bulk propagators of the orthonormal scalar fields fa satisfy the usual Klein-Gordon equations in AdSd+i 
space, as discussed in the previous sections. Explicitly, we have, in Euclidean signature and in the n — > limit, 

x (55) 

+ J- l/ (Ax)J l> (Ay)-2coa(vir)J. v (Ax)J- v (Ay)U, 

and 



(&#(?, v ,h) ' A<IA 



A 2 



' ' " + " ! ; (J„(Ax)J- v {Ay) + J_, y (Ax) J v (Ay) - 2cos(vir)J^(Ax)J- l/ (Ay)) 



r(l-^)(2A) 2 " 



(56) 

The conductance in Lorentzian signature, related to (ji(-p)jjip)) and (ji(—p))(jj(p)), are thus given by precisely 
the same 1-loop calculation as detailed in the previous section, except the scalar bulk-to-bulk propagators are replaced 
by G" and G 1 ^ 1 respectively. 



Given the notable complication in the form of the propagators ( 54 ), one would like to know if the divergence occur- 
ring in the collinear limit of the momenta, discussed extensively in the previous section, could actually acquire some 
n or boundary coupling / dependence, rendering our procedure in extracting the transport properties pathological. 
Fortunately, one can show that just as before, the divergent terms that are not extinguished by the phase space 
volume factor are again n and / independent. 

To see that, let us rewrite for example the basis rotation: 



-7=4>n ^ \ 0n-l = cos£0„ - sin£<£ n _i, (57) 

Jn V n 
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where we remind our readers that the tilde fields are those that diagonalize the boundary conditions. We deliberately 
rewrite the coefficients in terms of sine and cosine to make it explicit that the squares of these coefficients add up to 
one. The corresponding propagator is 



G" 



cos 2 £G nn + sin 2 £,G 



(«-!)(«-!) 



and to be explicit, 



("'" - I dk^ k2 Jv,nfJv,nf> G {n 1)( " ^ - j dA A 2 + k 2 J "fiJu,0, 



(58) 



(59) 



where J„,/ are defined in Ety. Plugging the propagator G nn into the loop including the external photon propagator 
then leads to the products 



x d/2 y d/2 K±_ 1 (-ipy)K£_ 1 (ipx)]J[cos 2 €J Utn f(A i x)J Ujn f(A i y)+sm 2 ^^ 



i=i 



2 r 



x (d-3)/2 y (d-3)/2 exp ( ipx ) exp(-ipy) J] 



7f 7f 

cos 2 ^(cos(A,x - 9 Vtn f - -))(cos(Ajj/ - 9 u , n f - j)) 



+ sin 2 ^(cos(A l a; - 6>„,o - -))(cos(A{j/ - 0„,o - -)) 



(60) 



where we have omitted some overall numerical factors and powers of Aj's which are common to all scalar loops 
independently of boundary perturbations. 

Concentrating on the collinear limit where one actually encounters a divergence i.e. Ai + A2 is close to p, the 
dominant contribution from the above expressions, after doing the radial integral over x and y is given similar to the 
discussion leading to (32), by 



(P - Al - A; 



^(cos^e-^^-^e^^t-^+sm^e-^ 



2e„,o-i) e »(2e„,o-f) 



+2 cos 2 £ sin 2 £e~* 



f+6~,o-f)pi(8„ int +9 u , -f) 



(p-Ai 



A 2 ) 1 - d (cos 2 £ + sin 2 2 = (p - Ai - A 2 ) 1 - d . 



(61) 



The dependence on n and / drops out in the leading divergence. One could repeat the exercise for the next leading 
order divergence to see that again it is n and / independent. Since the expansion in powers of n is completely regular 
in the propagators, the divergences remain n, f independent even as we expand first in n before doing the integral. 
We will in fact work out in detail in the following section an analytic example where the result after doing the radial 
integrals is particularly simple, and explicitly in this example the divergence can be seen to be n, f independent as 
we have demonstrated via the asymptotic expansion method here. 

We extract the dissipative part of the conductivity from the correlation functions via the Kubo formula 



1 



CT (Po) = — Im[(j l (-p )Ji(Po))retarded] 

Po 



(62) 



We include in figure (13) a numerical plot of the conductivity against coupling / at d = 4, v = |. The resultant 
conductivity is regulated, for convenience, by subtracting off a corresponding scalar loop satisfying Dirichlct boundary 
conditions in a simple un-replicated AdS. One can see that the conductivity decreases as disorder is turned on. This 
point shall be discussed again in our analytic example in the next section. 

We point out that if the random disorder has a non-trivial power function in momentum space, which corresponds 
to spatial corelation of these impurities, one can replace / by an appropriate function f(p) and the results are 
straightforwardly generalized by the same procedure. 



B. An analytic example in d = 2 + 1 



At d = 2 + 1 and v = ~, the result is particularly simple. The Euclidean photon propagator given by (16) for d = 3 
reduces to 



M*>p) = j u(p) 



(pzy/^Kiipz) 

(pe) a Ki (pe) 



9»J^(p)=0, K h (pz) 



(63) 
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FIG. 3. Dissipative part of the 1-loop contribution to the conductivity against coupling / at d — 4, v — i. External momentum 
being an overall scale is set to one. 



Similarly, scalar propagators proportional to products of J ± i (Ax) become simply products of y/2/irAx sin Aa; and 
■y/2/7rAa; cos Ax. More explicitly, (55) becomes in this case 



G u (x,y,k) = 
F u (x,y,A) = 



dA- 



A 



A 2 + k 2 
2cos(a;A) cos(yA) 
irA^/xy 



F u (x,y,A), 



I 



2itA 2 Jxy 



(cos(xA) sin(yA) + sin(ccA) cos(yA)) 



(64) 



Plugging into the loop integral and denoting the radial coordinates of the vertices x, y, in Lorentzian signature, we 
have 



(65) 



/disorder = / dxdyKi(ipx)Ki(-ipy)F nl (x,y,Ai)F nt (x,y,A 2 ) 

= 8(-l + Aj + A 2 ) 2 - f 2 (Aj + (-1 + A 2 2 f - 2A|(1 + Af)) 
VttAi A 2 (Af + (-1 + A2)2 - 2Af(l + Al))2 

where again we used the dimcnsionless quantities 



A, 



A, 
P 



f = 



!_ 
p' 



(66) 



It is important to note here that there is again a divergence not cured by the phase space factor in the collinear 
limit proportional to (p — Ai — A2) 1_d = (p — Ai — A2) -2 . However, it is completely independent of the boundary 
perturbation coupling /. Given the explicit analytic result of the integral, we shall subtract off the divergence directly 
here. 

Now we can compute the dissipative part of the conductivity in the presence of the disorder and find 

1 



o-(po)i-ioo P = ^-hn[{ji(-po)jj{po))i 
Po 



lim eS. 



(2tt) 3 p8 



loopJ 

;i i /•! 



o lJ 4(4Tr) 3 \Ki(i Po eW 



dAi 



(i 



l-Ai 



dAiAiA 2j ff(Ai,A 2 )J dis 



rdcr 



lim eSi 



2(£) 2 + W"^+log(8) 



y P0' 



Z2p Q {^Y\K h {ip e) 



(67) 



The /-independent part is simply the 1-loop correction to the conformal result. The precise value is unimportant, 
since it depends on how we subtracted the divergence. The /-dependent contribution however is interesting. First 
of all, it is negative definite, which means it always reduces the conductivity. Secondly, while the loop correction 
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is only of order N°, it is clear that in the low frequency limit (po/f 2 ) 2 <C iV~ 2 and this loop correction would 
eventually become more important than the tree-level contribution, which is of order iV 2 . This is not surprising since 
the boundary perturbation is deliberately chosen to be relevant in the infrared. This means that to extract result in 
the deep infrared limit one probably needs to re-sum the loop contributions. This could be the explanation as to how 
the delta-function at p$ — is corrected even though naively the boundary perturbation only begins to contribute at 
1-loop level and appears extremely weak. We would like to pursue this in more detail in future work. 

Time independent disorder 

When the disorder is time independent, only the component of the scalar 4>(k, z) with k = (0, k{) acquires a nontrivial 
boundary condition (T6J) . i.e. the disorder coupling / is frequency dependent and is more precisely given by fS Po fi. 
The disorder therefore only contributes non-trivially if at least one of the two internal bulk-to-bulk propagators is 
evaluated at exactly zero energy. By the Cutskovsky rules, the imaginary part of the loop is only non-zero when 
Ai ( 2 can satisfy the "on-shell" conditions that restrict them to be real and non-negative. All these conditions can be 
satisfied only at one point of the phase space where ft, = Ai = 0. In other words, by putting the KK-scalar fields 
on-shell and requiring that one of them has zero energy, its corresponding spatial momentum and effective mass has to 
be zero. Note however that the boundary disorder coupling / always come in the dimensionless combination (A) 2y //. 
Therefore as soon as we set A = such that the momentum dependent coupling / contributes, its contribution is 
immediately killed by the {K) 2u factor, for any positive v. The time-independent disorder therefore doesn't seem to 
have any effect on the dissipative part of the conductivity, even though the operator is relevant. The only possibility 
of a non-trivial contribution seems to be the marginal limit where v — > and Aq — > |. The computation should 
follow the same logic as presented in the previous sections, although the radial wave-function of the scalar field would 
involve Jo (Ax) and the boundary expansion (pi) has to be modified to include log-terms, in which case the correct 
linear combinations of Bcssel functions satisfying the relevant boundary conditions have to be computed separately. 
We will leave these interesting possibilities for future work. 

V. CONCLUSION AND DISCUSSION 

We have studied in the AdS space the 1-loop correction to the boundary-boundary2-point function of a 17(1) 
gauge field coupled to a complex scalar. The scalar is subjected to quadratic boundary perturbations and satisfies 
mixed boundary condition. These boundary perturbations are related to double-trace perturbations in the dual CFT. 
While we concentrated our effort toward studying the imaginary part of the diagrams, which is of particular interest 
physically and can be readily extracted via a simple generalization of Cutkosky rules [33J , we found that in Poincare 
coordinates these loops suffer from extra divergences shared by both the real and imaginary parts of the diagram, 
arising when the vertices of the diagrams are pushed toward the AdS horizon simultaneously and the geodesic distance 
between them shrinks to zero. When the loop integral is reduced to a integral of both vertex positions along the radial 
direction which gives rise to a product of two double variable hypergeometric functions (the Apell function F4), the 
divergence manifests itself as singularities along the boundary in the phase space where the momenta of the photon 
and the loop scalar fields become collinear. 

Fortunately, we manage to derive an explicit representation for these divergences and discover that they are inde- 
pendent of the boundary couplings, or in other words, the mixed boundary condition for the complex scalar. While 
F 4 is known only as a power series, and the singularities that concern us occur precisely along the boundary of its 
convergent domain, we were able to extract the exact properties of its singularities nevertheless using simply the 
asymptotic expansion of the radial wave-functions of the propagator. It is found that the order of divergences are 
finite and is related to the spacetime dimension in a simple algebraic way. The coefficients are easily computed as 
well. Most curiously, the leading divergence turns out to be universal, depending only on the spacetime dimension, 
whereas the first sub-leading divergence depends upon the spacetime dimension and the bulk mass of the scalar field 
only. The method presented can be applied to analyze the singularity structures of other functions that have a similar 
integral form, such as similar vertices in other curved spacetime background. 

These results are readily applicable to certain physical problems. In particular, we applied AdS/CFT correspondence 
to condensed matter systems and studied the effects of random disorder on conductivity, using the replica trick |10) . 
The effect of the disorder in this setup is captured by double-trace boundary perturbations, corresponding precisely 
to mixed boundary condition on the AdS side for the complex scalar. The dissipative part of the current-current 
correlation in the CFT is given by the imaginary part of the photon boundary-boundary correlator in the AdS space, 
which begins to exhibit the effect of the disorder from the photon coupling to the particular charged scalar at 1-loop. 
Since the singularities of the relevant loop integral are found to be universal, we are able to subtract the results of 
the conductivities obtained at various disorder coupling by the corresponding conformal result, and obtain a finite 
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difference. We found that in general the disorder reduces the conductivity. Also, since the disorder we consider 
is relevant in the infrared, we find that generally in the extreme low frequency limit, the loop correction begins 
to overwhelm the planar contribution, suggesting the need for re-summation, but allows for the possibility of the 
removal of the delta- function in the conductivity known to exist at zero frequency due to momentum conservation^ 
even though the effect of the disorder appears at first sight to be loop-suppressed. 

Let us mention a couple of words on some other interesting observations and conjectures. We found that the 
same divergence in the limit as loop vertices approach the horizon is quite general in loop integrals in AdS space, 
particularly when fields of higher spins, such as photons and gravitons are involved, although they most likely require 
more intriguing regularizations. However, for loops involving only scalar <p n vertices, the absence of divergence is 
correlated to renormalizability of the theory concerned. It is also interesting that the photon-fermion-fermion vertex 
displays a singularity near the horizon with precisely the same power law behavior as the photon-scalar-scalar vertex. 
It would be interesting to understand the precise physics of these divergences. In the Poincare patch where a mass 
gap is absent, we have, correspondingly in the bulk a continuum spectrum of "Kaluza-Klein" mass A that can be 
arbitrarily small, which might be intimately related to the collinear divergence in the decay amplitude of the photon 
as we have observed. 

Yet, it is suggested in |37| that in real-time computations one should build wave-functions such that they vanish at 
the horizon. This would mean that waves could only propagate as wave-packets rather than as momentum eigenstates 
along isometry directions. It is probable that the divergence can be cured by building a wave packet that vanishes in 
the horizon limit but asymptotes to a plane wave near the boundary. However, since the regular part of the diagram 
in the collinear limit is heavily suppressed by the phase-space volume, we expect that our results should not be too 
sensitive to the near horizon behavior of the wave-functions. We leave these important and rigorous endeavor for 
future work. 
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Appendix A: Contour integral of the scalar bulk-to-bulk propagator 

In the Euclidean signature, the scalar bulk-to-bulk propagator is given by: 



+ OG 



AdA [J v {Ax) + fA^J_ v {Ax)][J u {Ay) + fA^J_ v (Ay)} 



A 2 + k 2 



l + / 2 A 4 " + 2/A 2y cosz/7r 



(Al) 



We denoted / = (2A) 2l T(l - v)/[fT(l + v)] for brevity. 

We can evaluate this propagator by replacing it by a contour integral. Let's assume x > y for the moment. Notice 
that 



J u {Ax) + JA 2v J_ u {Ax) = 



(A2) 



Til 



and 



1 + / 2 A 4,y + 2fA 2lJ cos vtx = (/A 2 " + e v ™){fA 2v + e~ vlTl ) , 
so the original integral is also given by 



+oo 



AdA 



7ri(A 2 + k 2 ) 



K v {—ihx) 



K v (iAx) 



f A 2v p — wni/2 _i_ pVKijI f \2is pisiri / 2 _i_ p — VTzi/2 



(y terms) 



(A3) 



(A4) 



We define 



and 
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+°° AdA KJ-ikx) , 

(y terms), (A5) 



7Ti(A 2 + fc 2 ) jh2v e -vni/2 _|_ gi/Tri/2 

/•+ 00 AdA KJihx) , 

I-=- I —TlTo T^-~ ^7 --y terms. (A6) 

J 7Tz(A 2 + fc 2 ) f^2v e v7Ti/2 _|_ e -VTTi/2 yJ ' y > 

Let us also denote C+ the closed contour in the complex plane that consists of the entire real axis and the infinitely 
large semicircle in the upper half plane, and C_ the opposite contour that consists of the real axis and the infinitely 
large semicircle in the lower half plane. 
It's readily verified that if one substitutes 

A -)■ Ae" (A7) 

one finds 

1+ -*■ -I_ . (A8) 

Similarly, upon substituting A by Ae -7 ™, one finds 2_ becomes — J+. The two phase choices just mentioned, however, 



are mutually exclusive and therefore the full integrand (A2) does not enjoy any symmetry property as A — > —A. 



Let's evaluate the following contour integral, making use of equation (A7) and (A8| 

zdz KJ—izx) . 

,-a tttv -= , r ' (y terms) 

c+ lTl(z 2 + fc 2 ) f z 2v e -viri/2 + e »7ri/2 w 

AdA KJ-ikx) 

(y terms) 



fl+ 7Tl(A 2 + fc 2 ) f]\2v e -vni/2 + gViri/2 



AdA ^(-iAx) ( A9) 

K _ 7Tl(A 2 + fc 2 ) f\2v e -vivi/2 _|_ e viri/2 W > 

AdA KJ-ikx) , . 

(y terms) 






C+ _ J? 7Ti(A 2 + fc 2 ) j^2v e -VTxil2 + gi/TTj/2 



which gives rise precisely to the full integral | All. Here R + , R and R denote the positive, negative, and the full real 



axis respectively. The choice of the contour C+ is dictated by equation ( A7 1 as the phase of A rotates smoothly from 
to 7r. This contour integral can be evaluated, on the other hand, by picking up the residue at the pole A = fee 7 ™/ 2 . 
Similarly, we can evaluate the contour integral along C_: 



/ zdz Kjizx) . . 

J c _ 1Ti(z 2 + fc 2 ) f z 2v e vni/2 +e -vTri/2 ^ TU1) 

(y terms) 



JC_ Ki(z 2 + fc 2 ) f Z 2u e UTTl/2 _|_ g-I/TTJ/2 

AdA K u (ikx) 



R+ 7Tl(A 2 + fc 2 ) fJy2u e um/2 _|_ g-uiri/ 

AdA Kjikx) , , (A10) 

• (v terms) 

R _ 7Tl(A 2 + fc 2 ) f\2v e viTi/2 _|_ g-t-TTi/2 W ' 

AdA KJihx) , , 

(y terms) 



>C+-R Ki(k 2 + fc 2 ) jl\2u e PTTi/2 _|_ e iy7ri/ 

=1+ + 1- . 



This contour integral picks the residue at the pole A = fce -7 ™/ 2 and leads to an identical result. 

If x < y, both contour integrals given above diverge, in which case we should just exchange the role of x and y and 
carry out the above analysis similarly. 



In summary, we find the propagator ( Al ) equals 

2KJ[kx) \l u (ky) + fk 2 »I_J[ky) Q(x - y) 
U+I- = — 



1 + /fc 2 



+ x*+y. (All) 
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Appendix B: Some more details on the singularities of I v j x ,f 2 



We want to study the divergences of the integrals of the form [28] 



I 







z 2 K d _ 1 (iz) J ^{h-iz) J v {h.2z)Az 



T\l + \(ii + v)]T\\{d + n + v)}A%KZ 1 i , s a2 j2i 



(Bl) 



when 1 — | Ai| — | A.2I — > 0. Within the region |Ai| + | A.2 1 < 1, the Appell hypergeometric function F4 has a series 
expansion given by 



F i [a,b;c,d;k\,kl] = J2 



{a)m(b) r 



m!n!(c) m (cZ). 



A 2m a 2n 



where the notation ( • ) m is defined by 



(a), 



_ F (a + m) 

= r(o) ' 



(B2) 



(B3) 



This expansion fails to converge whenever 1 — |Ai| — | A.2 1 — > 0. One can, however, by choosing the integral contour of 
z carefully, extend the definition of / to the full (Ai, A2) plane which only becomes singular when 1 ± Ai ± A2 — > 0. 
To isolate these singularities of /, we examine the asymptotic behavior of its integrand as z — > 00 where the Bcssel 
functions have an asymptotic representation given by 



Kd_ x {iz) 



OO 



E 



M*) 



1 

2nz 



i(z-i 



E 



Here we defined notation 



(2z) n 



(v,ri) 



n=0 
2 4 ) 



(2iz) n 



E 



i n (2z)™ 



n!r(i + t/-n) 



(B4) 
(-7T < argz < 7r). (B5) 

(B6) 



and it is useful to note that (^,0) = 1 and (y, n) = (— z/, n). These expansions are only valid in the z-plane when its 
phase angle is between — -k and -k. Using these formulae, we can easily derive an asymptotic expansion for the full 
integrand of / in the limit z — > +00 and formally evaluate the integral I by integrating the resultant expansion term 
by term: 



1 



oo n n—k 



'■* (f-l,n- fc -Z) (/J, *;)(!/, Z) 
2f^A 1 Ai^ ^ ^y Z (2zY 



EEE 



^A^A 2 



. f e -i(e„+e v + s i-) i 2(k+l)-n e -ill-(A 1 +A 2 )]z + e _i(fl M _e„+f ) i 2k-n e -i[l-(A 1 -A 2 )]> 
+ e -i(e v -e„ + f) i 2l-n e -i[l-(A 2 -A 1 )]x +e i(8„+8 v + ^) r n e -i[l + (A 1 +A 2 )]x\ 



y y y (|-1."~ fc - (m, *0(^ 0W - n) 



(B7) 



2 n Aj=A 2 



n+ J 



2V7rAiA 2n=0fc=0I=0 

. (e-^fl^+^+^r^) (_)*+* [1 - (A! + A 2 )]" +i ^ + e -*(^-»-+¥) (_)*[i - (A x - A 2 )] 

+e- i ^-^+T)(_)'[i _ (A 2 - AO]^ 1 ^ + eWM+^+T^p. + (A : + A 2 )]" +i ^) , 
where we have defined the phase angle 



/(7T 



(B8) 
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As explained in details in section III C the "~" sign above must be understood as indicating that both sides equal 
apart from unknown finite functions. Therefore, terms that are not singular on the right hand side of the "~" sign 
are meaningless and should be discarded. In the domain we are interested in, < A 12 < 1, and the poles of / are 
given by 



-2i8„ r fd-1* 



i x i ,. \ i . v . i - e "" " r (V ) J 1 _ . 

v '°° y ' V 2iV^K^K 2 \[l-(A 1+ A 2 )]^ 



d-1 (i/,l) / 1 1 
4 d-3 l Ai A 2 



1 

[1-(A 1+ A 2 )] ! 



(B9) 



+ O ( [1 - (A x + A 2 )]" (d " 5)/2 J + finite things . 
More generally if < Aj. 2 < +00, I may contain other poles. Say, if 1 ± (Ai — A 2 ) — > 0, we find in the same way that 



Iv,oo(i,Ai, A 2 ) ~ - 



w: 



2 v /2 7 rA 1 A 2 1 [1 _ (Ai - A 2 )]^ 



d-1 (i/,l) / 1 



d-3 V Ai A 2 



1 



[l-(Ax-A 2 )]^ 



(BIO) 



+ o([l- (Ai + A 2 )]-( d - 5 )/ 2 ) + Ai <-> A 2 + finite things . 



Since the asymptotic expansion (B4| and (B5)are not valid along the negative axis of z, this method can not be used 



to analyze the pole structures if either Ai i2 becomes negative. In particular, the pole of the form of 1/(1 + Ai + A 2 ) a 
can not be properly derived as above, but they can be inferred easily by the symmetry properties of / when A — > —A. 
If, instead of the standard Bessel functions J±„, we wish to use J v j as defined by Ety in the integral formulae for 
/ as we should if the bulk scalar (f> satisfies a mixed boundary condition, we must replace the asymptotic expansion 
of J±„ by that of J u j- It's easily verified that 



J yj( z ) 



1 

2ttz 



P «(z-0i/,,f,A-f) 



71 = y ' 



E 



,,-i(z-e„ i/iA -j) 



E- 



(l*,n) 



L^, j»(2z)" 

n=0 v ; 



(Bll) 



which differs from those of J± u by merely a less trivial phase angle that we have denoted as 6vj,A above and is defined 
implicitly through the following equation: 



1 



tan 1 



W,A 



= tan 



a 2 " r(i-i/) 

, , a^ r(i-^) 
"*" f r(i+i/) 



tan#„ 



(B12) 



Obviously 9 v j interpolates between 9 V and 9- v = —9 V as / varies from 00 to 0. 

Therefore, it is straightforward to generalize the results given above and find, for example, when 1 — (Ai + A 2 ) — > 0: 



^.A./afrA^Aa) 



i(0u,/ I ,A 1 +0„,/ a ,A 2 )r(^i) 

2iV27rAiA 2 

d-l (u, 1) / 1 



[l-CAi+Aa)] 2 ? 



d-3 V Ai A ; 



1 



(B13) 



[1-(A 1 + A 2 )] ; 



O ([1 - (A! + A 2 )]- (d ~ 5)/2 ) + finite things 



which is essentially identical to that of J !/]OC except that 9 v j,h depends both on / and A now [38]. This is, of course, 
necessary if any nontrivial results are to be expected. 
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